A measurement-induced optical Kerr interaction 
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We present a method for implementing a weak optical Kerr interaction (single-mode Kerr Hamil- 
tonian) in a measurement-based fashion using the common set of universal elementary interactions 
for continuous-variable quantum computation. Our scheme is a conceptually distinct alternative 
to the use of naturally occurring, weak Kerr nonlinearities or specially designed nonlinear media. 
Instead, we propose to exploit suitable offline prepared quartic ancilla states together with beam 
splitters, squeezers, and homodyne detectors. For perfect ancilla states and ideal operations, our 
decompositions for obtaining the measurement-based Kerr Hamiltonian lead to a realization with 
near-unit fidelity. Nonetheless, even by using only approximate ancilla states in the form of su- 
perposition states of up to four photons, high fidelities are still attainable. Our scheme requires 
four elementary operations and its deterministic implementation corresponds to about ten ancilla- 
based gate teleportations. We test our measurement-based Kerr interaction against an ideal Kerr 
Hamiltonian by applying them both to weak coherent states and single-photon superposition states. 

PACS numbers: 



I. INTRODUCTION 



Optical nonlinearities, such as the Kerr effect, have 
great importance in both fundamental and technological 
sciences. The conventional way of obtaining an optical 
nonlincarity is through an intensity-dependant interac- 
tion between light and matter For instance, a self- 
Kerr interaction naturally occurs in an optical glass fiber, 
inducing an intensity-dependant phase shift upon a prop- 
agating light pulse. However, such natural Kerr nonlin- 
earities are extremely weak and a possible accumulation 
of the nonlinear interaction through a sufficiently long in- 
teraction time, e.g. via a long-lasting pulse propagation 
over a long fiber segment, is impossible due to the domi- 
nating effect of accumulating linear losses (however, see, 
for example [|J]). Thus, research on optical nonlinearities 
is typically based on having an exotic medium with spe- 
cial optical properties in order to realize a more robust, 
though still usually fairly weak nonlinear interaction be- 
tween the medium and the light field. This research field 
is very active and produces various promising results [3|- 
5] . Most recently, rather strong Kerr- type couplings were 
achieved in electro-optical circuit QED systems Q. 

There is also the more recent approach of 
measurement-induced nonlinearities. The most promi- 
nent example for this is an efficient linear-optics 
quantum computation scheme that relies on entan- 
gled multi-photon ancilla states (in order to become 
near-deterministic) and photon counting 0, H[. This 
approach works on the single-photon level and uses 
discrete-variable qubit-type encoding. However, there 
are also measurement-based proposals to implement 
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a so-called continuous-variable (CV) cubic phase gate 
(the hardest gate to implement from the CV universal 
gate set), which, in principle, can b e ap plied upon an 
arbitrary optical single-mode state 0, Il0j |. A very recent 
result is the experimental demonstration of such a (weak) 
cubic phase gate 11] . This particular demonstration is 
the key to engineer an arbitrary nonlinear interaction, 
since any other elementary CV gate besides the cubic 
phase gate is a Gaussian operation and hence fairly easy 
to realize [lU, [l3| . However, current theoretical work on 
using measurement-induced, elementary nonlinear gates 
to obtain more useful nonlinear Hamiltonians, such as 
the Kerr effect, has been still insufficient with regards to 
an actual experimental implementation. 

In this work, we present an alternative method to the 
medium-based Kerr nonlinearities and demonstrate ob- 
taining the Kerr nonlinearity in a measurement-based 
fashion within a few measurement steps. Our imple- 
mentation of the Kerr gate is independent of the input 
state and works, in principle, for arbitrary optical states. 
Interestingly, in contrast to the medium-based schemes, 
our approach works best on the level of individual pho- 
tons. Our scheme could be more generally referred to 
as measurement-induced Hamiltonian engineering, how- 
ever, here we shall only focus on the important example 
of the Kerr interaction. 

The main theoretical background consists of the re- 
cently developed nonlinear-gate decomposition idea [l3j 
combined with an approximate implementation of the cu- 
bic phase gate 10]. However, besides combining them, 
we also extend these two ideas, presenting a rather new 
approach in order to avoid any redundant steps for both 
decomposing the Kerr interaction and preparing the an- 
cilla states. As a result, we will be able to realize the 
Kerr gate using a very small number of steps (that is 
between 4 and 10 steps, as opposed to the tens of steps 
needed in the original scheme of Ref. [l3[ ). For implc- 
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meriting the elementary interactions, which are slightly 
more general than the common universal CV gates, we 
apply the measurement-based (cluster-based) model pro- 
posed in Refs. [3 EH- For this purpose, we first show 
that four elementary interactions are enough for realiz- 
ing the full sclf(singlc-modc)-Kcrr interaction with an in- 
teraction strength around the same magnitude as the 
medium-based Kerr nonlinearities. However, in order 
to implement these four interactions in a measurement- 
based and deterministic fashion using elementary CV 
cluster teleportations, six additional correction opera- 
tions are needed. By incorporating these corrections into 
the measurement-based protocol, in total ten teleporta- 
tions will be required. Since the elementary operations 
are of cubic and quartic order in the mode operators 
and hence correspond to non- Gaussia n op erations, in the 
original CV cluster-based model [3, |T^, the correction 
operations would correspond to quadratic squeezing and 
further non-Gaussian operations, which must be taken 
into account in the choice of the measurement bases at 
each step of the cluster computation. In our ancilla-based 
model, however, the measurements are always homodyne 
detections, and for the at most cubic correction opera- 
tions, suitable cubic ancilla states have to be prepared 
depending on the previous measurement steps. When the 
approximate, heralded single-photon-based ancilla states 
are employed [lOj for this purpose, it would mean that 
the quantum information to be processed in the clus- 
ter has to be temporarily stored in a quantum memory 
13]. Without such a memory our proposal would be 
probabilistic. However, if one has access to an efficient 
quantum memory, our weak-Kerr scheme could be also 
concatenated in order to implement a strong Kerr in- 
teraction. These concatenations would accumulate the 
errors from each weak-Kerr scheme though, and thus, for 
a strong-Kerr scheme to attain still high fidelities, the 
gate decompositions of the weak-Kerr schemes have to 
be refined through higher-order approximations, as we 
will show. 

This work is organized as follows: in Sec. [Til w e de- 
scribe the measurement- based model of Refs. |l4l. fl5| . 
Then, in Sec. IIII1 we present a new approach to decom- 
posing a nonlinear interaction and apply it to the Kerr 
interaction. After this, in Sec. lIVi we will describe an ex- 
perimental proposal for the necessary ancilla states and 
provide an explicit recipe for the implementation of the 
Kerr interaction. In Sec.fVl we test our scheme by apply- 
ing the resulting decompositions to different input states 
and comparing the outgoing states with those obtainable 
via the exact evolutions. Then, in Sec. I VII we present 
a guideline to obtain higher-order decompositions and, 
finally, in Sec. I VIII we briefly discuss the possibility of a 
strong Kerr interaction. 

Throughout, we use capital letters and hats for opera- 
tors and small letters for scalars and functions. In equa- 
tions, the letter i is used only as the square root of -1. 
The notation [A, B] is used for the commutation of the 
operators A and B. We use the convention h — 1/2, i.e., 



the fundamental commutation relation is [X, P] = i/2 
with X = (at +a)/2 and P = i{a) - o)/2. 



II. MEASUREMENT-BASED 
IMPLEMENTATION 

For the experimental implementation of the required 
operators, we use the teleportation-based model of 
Refs. [l3, Us} . According to this model, it is possible 
to teleport elementary interaction Hamiltonians (gates) 
that are diagonal in X, including rotated versions of 
them, onto an arbitrary input state using Gaussian cou- 
pling, homodyne detections, and an appropriate set of 
ancilla states. 

Assume that we have the (single-mode) state \ip) that 
we want to transform, and the operator we intend to 
apply is A(X), which is an operator diagonal in X. We 



aim to obtain the state 



' final 



A(X) where 



ip(x) \x) dx , 
A(X) \x) = a(x) \x) , 

final = J Oi{x)il{x)\x) dx . (1) 

In order to obtain the state W) fi na ii we couple the input 
state \ip) to an offline-prepared (single- mode) state \a) 
using the Gaussian coupling operator e 2lX ^ x ^ ; where the 
subscripts refer to the two modes. This coupling can be 
exactly implemented by beam splitters and squeezing op- 
erations [161 ] . Then we make a P homodyne detection on 
the input mode (mode 1), which yields the state \^) output 
in mode 2 with 



| a) = / a(x) \x) dx , 
W)o. P . = 7^/(/ e 2iy ^-^(y)dy)a(x)\x)dx(2) 

where f3 is the measurement result of the homodyne de- 
tection. Equation ([2]) can be written as follows, 



W output = A (X)F j e 2if>X W 



(3) 



Here, the operator F is the Fourier transform operator 
defined as follows, 



F\x) = 



1 



2ixy 



\v) dy ■ 



(4) 



which can be used to transform the operators in X to 
operators in P and vice versa by unitary conjugation, 



Fe ltpm F f = e ( ' 1)mitxr ' 



(5) 



To obtain the desired state JT]) from the output state ([3]) , 
we need to apply further correction operations. For this 
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purpose, Eq. Q can be written as follows, 

{^output = A(X)F^e^ x FFU^) 

= A{X)e~ 2tf3P F^ \i/>) . (6) 
Then we implement a correction operator O correc , 

1^) final = O correc W output 

= A(X)F*\i,), (7) 

which coincides with the desired final state \ijj) final up 
to the inverse Fourier transform F\ and where 

O CO rrec =A(X)e 2ll3P A\X). (8) 

For example, when we want to implement the cubic phase 
gate, A(X) = e ltx , we need the following Gaussian cor- 
rection operator, 

e 2ipP+3itpX 2 /g-j 

For the quartic gate A(X) = e ltxi , we need a cubic cor- 
rection operator, 

2i/3P+4it/3X 3 _ -3it£ 4 p 2i/3P iitpX 3 -6it/3 2 X 2 4it/3 3 X 

(10) 

Note that the displacement operations, e ltx and e ttp , can 
be realized by just coupling a bright light beam with the 
signal state through a beam splitter. 

Now in order to get rid of the inverse Fourier transform 
in fi na i, one possible strategy would be to first apply 
a Fourier transform upon the input state, = F \ip), 
prior to running it through the above protocol. This 
eventually gives the desired output state, \ijj}f ina i = 

A(X)F* \fl>)' = A{X) |V) = m finaV Applying a Fourier 
transform in a measurement-based fashion is very easy: 
it just means teleporting the input state in a single, el- 
ementary step using one of the canonical ancilla states 
for approximate CV cluster computation, namely a P 
squeezed vacuum state 0, [l5[ . 

Another strategy is to leave the inverse Fourier trans- 
form in the final states of the protocol and modify the 
interaction sequence in accordance to this. For example, 
suppose we want to teleport the four operators A, B, C, D 
onto the initial state \ip) . We may then teleport the modi- 
fied operators A',B',C, D' onto the input and we obtain 
the final state, 

A'F^B'F^C'F^D'F^ \i>) 
= A'F^B'F^C'F^F^FD'F^ 
= A'F^B'F\F^f{FfC'{F^) 2 FD'F^ |V) (11) 
= A'{F^f{FfB'{F^f{FfC'{F^) 2 FD'F^ \ij>) 
= A'{FfB'{F^f{F) 2 C'{F^) 2 FD'F^ |^) , 

where 

D' = F^DF, 
C = (F^) 2 C(F) 2 , 
B' = (F^) 3 B(Ff, 
A' = A. 



While the operators A, B, C, D are assumed to be diag- 
onal in either X or P, the operators A',B',C',D' are 
straightforward to derive. 

All arguments used in this section for operators diag- 
onal in X also hold in a similar way for the operators 
diagonal in P. 



III. DECOMPOSING THE KERR 
INTERACTION 

The general problem in decomposing a given interac- 
tion is to implement the desired operator using a limited 
set of experimentally available operators. In our previ- 
ous work [IH , we presented a systematic and efficient 
recipe for this, however, in order to optimize these de- 
compositions even further with regards to an experimen- 
tal implementation, here we shall use a slightly different 
approach. We will avoid any specific operator approxima- 
tions such as splitting and (nested) commutator approx- 
imations, as employed in Ref. (l3| . and instead utilize a 
more brute force approach especially tailored for the Kerr 
interaction. This approach also has a greater flexibility 
for choosing the order of the operators, which turns out 
to be beneficial in reducing the number of experimental 
steps when additional correction operations are needed 
later on. 

Let us first define the Kerr interaction as 

e u(x 2 + P 2 -i) 2 = ^t{X 2 + P 2 f e ^t(X 2 +P 2 -\) 

_ e it(x 4 +p' l +x 2 p 2 +p 2 x 2 ) e -it(x 2 +p 2 -±) 

= e t^X* + l P* + ±[X\P 3 ]- t l s t) e -rt(X 2 +P 2 -i) 

oce ,(^ + iP*+S[^,p»]). (13) 

Here, at the beginning, compared to the standard self- 
Kerr interaction e 4 *^^ -1 ), with the number operator 
N = X 2 + P — i, we omitted the Gaussian phase- 
rotation operation e~ ltN ; and also in the last line we 
omitted that phase rotation and another global phase 
factor for simplicity. Gaussian operators are considered 
as relatively easy to implement and they have been al- 
ready achieved experimentally [l7j . Therefore, we define 
the Kerr interaction using the operators X 4 , P 4 , X 3 , 
and P 3 . In this form, and using decompositions into 
operators of this form, we can then directly apply the 
measurement-based model introduced in the preceding 
section, Sec. |TTJ Note that in our previous work [13j, us- 
ing our general and systematic decomposition scheme, we 
substituted X 4 and P 4 by nested commutators of X 3 , 
P 3 , X 2 , and P 2 , and we employed correspondingly a 
universal gate set with only quadratic and cubic gates. 
Here, we do allow for quartic interactions in the elemen- 
tary gate set, which nonetheless can be incorporated into 
the optical, ancilla-based implementation scheme. 

In order to avoid splitting and commutator approxi- 
mations, we aim to decompose the Kerr interaction using 
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concatenations of the following set of interactions, 



A. With postselection 



Taking the logarithm of arbitrary concatenations of the 
two elementary interactions in Eq. (|14|) leads to the fol- 
lowing elements, 



{tX 3 , tP 3 , tX 4 , tP 4 , t 2 [X 3 , P ], 
t 2 [X\P%t 2 [X\P\t 2 [X\P\...}. 



(15) 



We omit the higher-order nested commutations. Pro- 
vided that t <^ 1, it will be sufficient for us to only deal 
with these elements up to a negligible error. This means 
that we have eight conditions to satisfy. It is generally 
possible to reduce the number of operators in the decom- 
position by reducing the number of necessary conditions 
to satisfy. In order to do this we divide the operators X 3 
and P 3 by \fl in Eq. (TT4")) . Now, the concatenations will 
lead to the following elements. 



{t x ' 2 X 3 , t^ 2 P 3 , tx\ tP\ t[X 3 , P% 



(16) 



t 3/2 [X 3 ^ p4] > 1 3/2 [X A > pZ ] ^ [X \P%...}. 

Hence, it is sufficient to only satisfy the conditions for 

{t 1/2 X 3 , t 1/2 P 3 , tX 4 , tP\ t[X 3 , P 3 ]} , (17) 

neglecting any orders 0(i 3 / 2 ). Now consider the follow- 
ing concatenation, 



e ipit 1/2 P 3 +ip 2 tP 4 e i P3 t 1/2 X 3 +ip A tX l 
xe ip5t 1/2 P 3 +tp 6 tP i e ip 7 t 1 ^X 3 +ip s tX i 



Taking the logarithm of this concatenation yields 



(18) 



»(P3 + Pl)t 1/2 X 3 + i(px + p 5 )t^ 2 P 3 + i( Pi + P8 )tX 4 + 



i{P2 + Pe)tP A + ~ (pips - P3P5 + PxPi + P5P7) t[X 3 , P 31 



+ 0{t 3 ' 2 ) 



(19) 



We then solve the corresponding polynomial equations 
for the coefficients pt [l3| in order to obtain the logarithm 
of the desired Kerr operator in Eq. (fT3|) . One possible 
solution set for Eq. (fT9"l) is as follows, 



e tt^ 2 P 3 e lit^ 2 X 3 e -it 1 / 2 P 3 +itP i e - l it l ' 2 X 3 +itX i _ ^ 

These four operators, each diagonal either in X or P, 
can be implemented in a measurement-based fashion us- 
ing suitable ancilla states, as described in Sec. HU Before 
discussing how to possibly realize such ancilla states with 
quantum optics, let us address the question whether tele- 
porting the above four operations onto an arbitrary input 
state can be done deterministically. 



In order to avoid the necessary correction operations 
when implementing the non-Gaussian cubic and quartic 
gates in a measurement-based fashion (recall Sec.|ll|, we 
may simply postselect the measurement result j3 « 0, 
in which case each of the operators in Eq. (f20| corre- 
sponds to one teleportation step. Then we only need four 
non-Gaussian ancillae to decompose the Kerr interaction. 
Following Eqs. (jTTJ) and (TT2)) . the actual operators to be 
teleported onto the input state are given by 



a t 1 ' 2 x 3 + l tx' 1 , at 1 ' 2 P 3 +itp* 
e -.i 1 ' 2 p 3 + l iP 4 _^ e it 1 ' 2 p 3 +itp i 

e ^' 2 X 3 _^ e tit 1 "P a 
e it^ 2 P 3 _^ el t^ 2 P 3 



(21) 



B. Without postselection 

For a deterministic operation one needs to implement 
the corrections. In this case, we choose to realize each 
cubic and quartic interaction separately, 



j t l/2p3 *y 



1/2 . 



,1/2 



itX* 



(22) 



This concatenation is now constructed in order to be able 
to implement the correction operations for the quartic in- 
teractions at the same time together with the subsequent 
cubic interaction and its corresponding quadratic correc- 
tions. Recall that for the quartic operator we need cubic 
and quadratic corrections, while for the cubic operator 
we need quadratic corrections. Therefore, in total we 
need ten elementary teleportation steps. Displacements 
do not require any teleportation steps, thus we do not 
count them. 

One important subtlety, however, must be mentioned 
here, when the ancilla-based scheme of Scc.|TT]is to be em- 
ployed for realizing the ten teleportation steps. Since the 
correction operations depend on the homodyne measure- 
ment outcomes, the corresponding ancilla states cannot 
all be prepared beforehand, i.e., offline in a strict sense. 
For every gate that depends on a measurement outcome, 
the right ancilla state must be prepared effectively on- 
line. A possible remedy for this dilemma is the use of a 
quantum memory, in which the quantum information to 
be processed is temporarily stored until the conditional 
ancilla-state preparation (e.g. in the form of photonic 
superposition states, see later) has succeeded fl3j |. 

The scheme without postselection would be certainly 
most powerful, if even the non-Gaussian ancilla states 
could be prepared in a deterministic fashion whenever 
they are needed. In this case, ten deterministic tele- 
portations with ten deterministically prepared ancillae 
would give a fully deterministic Kerr gate. Compared to 
ten deterministic teleportations with six probabilistically 
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prepared ancilla states (corresponding to the six opera- 
tions in Eq. (j2"2")l while the extra four quadratic correc- 
tion operations rely upon unconditional squeezed-state 
resources), the postselected scheme with four probabilis- 
tic teleportations and four probabilistically prepared an- 
cillae may well be more efficient; however, the necessarily 
finite postselection window would also add extra errors 
in addition to those caused by imperfect ancilla states. 



IV. APPROXIMATING THE CUBIC AND 
QUARTIC INTERACTIONS 

In our measurement-based Kerr-interaction scheme, 
the ancilla states from Eq. @ that correspond to the 
necessary cubic and quartic operators are 



").,= / e ltx \.r 



itx 3 -\-it'x 



x) dx , 


(23a) 


x) dx , 


(23b) 


x) dx . 


(23c) 



Experimentally realizing the states in Eq. (|2"3")l is the most 
challenging part of our scheme. For thispurpose, we use 
the Taylor expansion approach of Ref. [lfj and extend it 
from cubic to quartic gates (interactions). The simplest 
approximations for cubic and quartic interactions and 
their corresponding ancilla states are as follows, 



After expanding the state (|26|) , we obtain the conditions 
for the state (|24a|) . namely 

ci + c 2 + c 3 = , 
3 

cic 2 + c 2 c 3 + cic 3 + - = , 
itc\C2Cz — 1 = 0. 

For example, for a cubic interaction with the amplitude 
10 -3 , the following set is one possible solution, 

ci 9.95i , 
c 2 -> -8.70356 -4.975i, 
c 3 -> 8.70356 -4.975i. 

The necessary procedures for the other operators can 
be derived in a similar fashion. 



V. EXAMPLES 

In this section, we apply the derived decompositions 
upon some test states and compare the results with those 
for the same states under an exact Kerr evolution. We 
first assume to have access to perfect ancillae. In other 
words, we assume that we have the elementary operators, 
such as e ltx +tt , available on demand with perfect 
fidelity. Later we shall relax this condition and assume to 
have access only to approximate versions of the ancillae, 
like in Eqs. (JUJ). 



\a) 1 m J (1 + itx 3 ) \x) dx , (24a) 
I a) 2 « J {I + itx i ) \x) dx, (24b) 
|a) 3 w J {I + itx 3 + it'x i ) \x) dx . (24c) 

Even though this first-order approximation is very lim- 
ited, in order to obtain a small-amplitude Kerr inter- 
action it is sufficient. Moreover, an approximate cubic 
interaction like in Eq. ()24a[) has already been observed 
experimentally (Tlj . 

The necessary ancilla states (f2~4"f can be created by 
applying photon substraction and displacement opera- 
tions upon a squeezed state similar to 10], [l8j]. Note that 
a suitable squeezed vacuum state S |0) can be approxi- 
mated by an eigenstate of the operator P with eigenvalue 
zero, \p = 0), such that 

D(-ci)a J D(c 1 )£>(-c 2 )a J D(c 2 )£>(-C3)a£>(c 3 )S|0) 

= (a + Cl )(a + c 2 )(a + c 3 )S|0) , (25) 

corresponds to 
{X + iP + Cl )(X + iP + c 2 )(X+iP + c 3 )\p = 0) . (26) 



A. Examples with ideal ancillae 

As an example for using the ideal ancillae, the accuracy 
of our Kerr decomposition is illustrated in Fig.Q] The ex- 
act Kerr evolution with amplitude 10~ 3 and the decom- 
posed evolution are applied to a coherent state which has 
an amplitude of ft = 1. In the position representation, 
this state is given by 

\P)= 0) ' J e-^ 2 \x)dx. (27) 

In order to be able to calculate the exact evolution, 
~ e ltN , one can use the well-known Fock expansion for 
the coherent state. Note that in Fig.[TJ both the real and 
the imaginary parts of the exact and the decomposed 
evolutions are overlapping, indicating the high accuracy 
of the decomposition (while there are no additional er- 
rors, in the measurement-based model corresponding to 
the use of an ideal ancilla-state set). 

The quality of the approximation depends on the am- 
plitude (strength) of the applied interaction as well as 
on the input state. For comparison, we define an er- 
ror in the decomposition through the inner product of 
the exactly evolved state with the output state of the 
same input state under the decomposed evolution: e = 
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Real pari of the exact evolution and the decomposed evolution Real part ot the exact evolution and the decomposed evolution 




FIG. 1: Kerr evolution of a coherent state, \/3 = 1). The 
position representation of the decomposed and the exact evo- 
lutions are shown. The blue line corresponds to the de- 
composed, the red dashed line to the exact evolution. The 
strength of the Kerr interaction here is 10 -3 . In the cor- 
responding measurement-based implementation, the ancillae 
are ideal. 



FIG. 2: Kerr evolution of a coherent state, \/3 — 1). The 
position representation of the decomposed and the exact evo- 
lutions are shown. The blue line corresponds to the de- 
composed, the red dashed line to the exact evolution. The 
strength of the Kerr interaction here is 10 _1 . In the cor- 
responding measurement-based implementation, the ancillae 
are ideal. 



1 — | (ipexactlipapprox.) |- We use the coherent state as a 
test input state and the error values for different config- 
urations are listed in Table |U As one can see from the 
table, an increase of the amplitude of the coherent state 
also increases the error. In general, our scheme works 
better for weak states, in the sense that the light fields 
contain only a small number of photons. Similarly, an in- 
crease of the amplitude (strength) of the Kerr interaction 
increases the error. Nonetheless, one can also implement 
a reliable, high-amplitude Kerr interaction, as we will 
discuss later. 



inter, amp./ coherent state amp. 


1 


5 


10" 3 


-^q-7.7594 


■^q-2.6641 


10 _a 


^q-4.9653 


■^q-0.7526 



TABLE I: Errors for different configurations. 

Note that there are also numerical errors for obtaining 
the decomposed evolutions, which contribute to the total 
errors. It was not possible for us to extract the numer- 
ical errors alone, and so we cannot present those errors 
which are solely due to the interaction decomposition. In 
this sense, one can infer that the actual accuracy of our 
decompositions is even better. 

As an example for a case with a significant error, we 
employed our decomposition for a Kerr amplitude of 
10 , still with a coherent state where (3 = 1. The results 
can be seen in Fig. [SJ The errors are now clearly visible. 

We also applied our decomposition upon a superposi- 
tion state of a single photon and a vacuum, 

-^(|0) + 11)), (28) 

for an interaction amplitude of 10 -3 . This state is sim- 
ilar to a very weak coherent state, oc |0) + /3 11), but as 



opposed to a coherent state it is highly non-Gaussian. In 
this case, the error was 1Q -7 - 9853 . 



B. Examples with imperfect ancillae 

When all the interactions in Eq. (|20|) are replaced by 
the operators (|24p , and again considering the previous ex- 
ample for the Kerr interaction with an interaction ampli- 
tude of 10 -3 applied to a coherent state with amplitude 
1, we obtain an evolution as shown in Fig. [3l Compared 
to the case of ideal ancillae, one can see that there is an 
increase of the error from io~ 7 - 7594 to io -2 - 7446 due to 
the approximation of the ancilla states. 



Real pan of the exact evolution and the decomposed evolution 




■0.5 1 1 1 1 1 1 1 1 1 1 

-10 -8 -6 -4 -2 2 4 6 8 10 

x iq- 3 Imaginary part of the exact evolution and the decomposed evolution 




-10 -8 -6 -4 -2 2 4 6 8 10 



FIG. 3: Kerr evolution of a coherent state, \/3 — 1). The 
position representation of the decomposed and the exact evo- 
lutions are shown. The blue line corresponds to the de- 
composed, the red dashed line to the exact evolution. The 
strength of the Kerr interaction here is 10 -3 . In the cor- 
responding measurement-based implementation, the ancillae 
here are only approximate. 
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VI. DERIVING HIGHER-ORDER 
DECOMPOSITIONS FOR THE KERR 
INTERACTION 



The last equation corresponds to the cross terms of the 
first and the second-order terms. One solution set for 
this equation system is the following, 



In order to be able to implement the Kerr interaction 
with a smaller error, we need higher-order decomposi- 
tions. In this section, we will present an efficient and 
direct approach to derive such higher-order approxima- 
tions. Let us start with some definitions. 



1. Using the method explained in Sec. IIII| one can 
check the following operator decomposition, 



e'2 



tX* AtP'+itP* 



tx-> 



itP gi^tX —i^tX 



^itX i +itP i + it 2 [X 3 ,P 3 ]+t 3 F 1 + t*F 2 + ...) 



(29) 



After the replacements X 3 — > ^ and P 3 -» 
this is a second-order decomposition for the Kerr 
interaction and we denote it as Q2(t). For our pur- 
pose, the specific forms of the operators F\, F 2 , ... 
are not important. 

2. The inverse of the decomposition in Eq. (|2U1) . 

can be obtained by reversing the decom- 
position and replacing t with —t, 

it* 4 _ 

~ (30) 



e -i±tX 4 +i±tX 3 e itP 3 e 



\tX 3 ^itP^-itP 3 - 



,(- 



itX" 



lt 2 [X 3 ,P 3 ]-t 3 F 1 -t' L F 2 + ...) 



This decomposition can be used for the Kerr inter- 
action with a negative amplitude. 

3. Thus, when we reverse the decomposition (f2l?|). 



Q 2 ev (t), we obtain the following operator, 



yrev 

ihtx 4 



i£tX 3 -UP 3 A%tX 3 JtP^+itP 3 „iitX a _ 
e (itX i +itP i ~%t 1 {X 3 ,P 3 ]+t 3 F 1 -t i F 2 + ...) 

4. Another useful decomposition is the following, 



(31) 



\tX* -UP* 



UP 3 e -i±tX 3 e itP 3 e ~i±tX i +i±tX 3 



D (-itX 4 -itP* + it 2 [X ? \P 3 ]-t 3 F 1 +t i F2 + ...) 



(32) 



Now in order to derive a third-order approximation, 
we will use a concatenation of Q2(t) and Q? e "(t). Similar 
approaches have been presented in Refs. [13l Il9|. Hence, 
an order condition for a third-order approximation can 
be derived with the help of the following concatenation, 



i{c 1 t)Q r 2 ev {c 2 t)Q 2 {c s t)Q r 2 ev {c i t). 



(33) 



Eliminating the third-order terms while keeping the first 
and second-order terms corresponds to the following or- 
der conditions, 



Cl + c 2 + c 3 + c 4 



4 + 4 + 4 + c \ 



1, 
0, 



(34a) 
(34b) 
(34c) 



c\c 2 + C X 4 + CJC3 - C2C3 - Ci4 - C2C3 + CJC 4 



2 1 2 1 2 
-C 2 C4 + C3C4 + C1C4 



C2C4 



C3C4 



. (34d) 



("2 
C3 

C4 




(35) 



and by making the replacements A 3 — > ^ and P 3 — > 

^t=, one can obtain a third-order decomposition for the 
Kerr interaction. Further higher-order decompositions 
can be derived in a similar fashion. 

When applying the above third-order decomposition 
upon a coherent state with amplitude j3 = 1 for a Kerr 
interaction with amplitude 10~ 3 , like in the example 
of Sec. IV Al the error in the final state decreases from 
2Q-7.7594 |- 20~ 9 - 9569 . An error reduction like this is 
crucial for concatenating weak Kerr gates (Kerr interac- 
tions with small amplitudes) sufficiently many times in 
order to obtain strong Kerr gates (Kerr interactions with 
large amplitudes), which we discuss in the next section. 



VII. HIGH-AMPLITUDE KERR INTERACTION 

We can obtain a high-amplitude Kerr interaction by 
applying a suitable decomposition many times. As an ex- 
ample, we apply our third-order decomposition for 1000 
times with an initial amplitude of 10~ 3 . This enables us 
to obtain a Kerr amplitude of I. We apply this whole de- 
composition upon a coherent state with amplitude /3 = 1, 
for which the results are shown in Fig. |4j The error for 
this particular decomposition is 10 -3 - 8252 . 

As another example, we apply the high-amplitude Kerr 
gate to simulate the effect of a nonlinear sign shift gate, 
which can be used to entangle photons [20j . More specif- 
ically, we apply the third-order approximation for the 
Kerr gate onto the following state: 



(|0> + |1> + |2», 



(36) 



which is supposed to transform to: 



(|0> 



|2». 



(37) 



For the initial strength of it x 10 3 , we apply the de- 
composition 500 times and find an error of 10 -3 2423 . 



8 



Real pari of the exact evolution and the decomposed evolution 




FIG. 4: Kerr evolution of a coherent state, \/3 — 1). The 
position representation of the decomposed and the exact evo- 
lutions are shown. The blue line corresponds to the decom- 
posed, the red dashed line to the exact evolution. A Kerr 
interaction with an amplitude of 10 -3 is sequentially applied 
for 1000 times in order to get a total Kerr interaction with an 
amplitude 1. In the corresponding measurement-based imple- 
mentation, the ancillae are again ideal here. 

VIII. SUMMARY AND CONCLUSION 

In this work, we presented an experimentally accessible 
method for implementing a weak Kerr interaction upon 
an arbitrary single-mode optical state in a measurement- 
based fashion. For this purpose, the Kerr interaction can 
be decomposed into a sequence of 4-10 elementary opera- 
tions. Each Hamiltonian for these elementary operations 
is either cubic or quartic in X and P such that the nec- 
essary ancilla states in the measurement-based realiza- 
tion are superposition states of up to four photons. The 
measurements are simple, efficient homodyne detections, 
while the signal and ancilla states are coupled through 
linear beam-splitter transformations. 

For the case that the homodyne measurement results 
are postselected around the origin in phase space, there 
are no correction operations needed, and so, four ele- 
mentary Hamiltonians suffice. In order to avoid such 
postselections, it is necessary to do the correction opera- 
tions, which can be of quadratic or even cubic order and 
will depend on some of the measurement outcomes. In- 
cluding these corrections, ten elementary interactions are 
needed, and correspondingly many ancilla states. Since 
the cubic ancillae can be optically prepared only in a 



heralded fashion, in this case, quantum memories would 
have to be employed. If the necessary ancilla states for 
the elementary operations including the corrections were 
all available on demand in a deterministic fashion, we 
would obtain a fully deterministic Kerr gate. 

We illustrated our Kerr decompositions with various 
examples, where test input states such as coherent states 
and single-photon states were subjected to an exact and 
an approximate, decomposed Kerr evolution. These ex- 
amples showed that our decompositions have a high ac- 
curacy, which mainly depends on the amplitudes of the 
input states as well as the amplitudes of the Kerr inter- 
action itself. The smaller these amplitudes, the higher 
the accuracies. This is why our scheme works best for 
states with fairly low (mean) photon numbers. 

In order to obtain a strong Kerr interaction, we also 
discussed how to concatenate our decompositions for 
the weak Kerr interactions and showed that for this 
purpose, higher-order decompositions are essential in 
order to avoid the accumulation of errors. In princi- 
ple, provided sufficiently many, suitable cubic ancilla 
states are available on demand, a fully deterministic, 
strong Kerr interaction would be implementablc in a 
measurement-based fashion, representing an alternative 
approach to medium-based nonlinear dynamics. 



Appendix: Exactly decomposing the quartic 
interaction 

It is also possible to exactly realize the operator e ltx 
by Gaussian operators and cubic operators, using the two 
relations, 

e |it?t 2 X=X 2 _ e -it 2 X| e it 1 X 1 P 2 e it 2 X% £ -2it 1 X 1 P 2 

xe i * 2X 2 3 e it i x i p 2 e -i*2^ 2 3 

and 

e it2-xf _ e itx 2 t p 2 e %tx\x 2 e -itx*x 3 e -itxfp 2 

where the subscripts refer to the corresponding modes. 
In a similar fashion, even arbitrary powers of X and P 
might be decomposed exactly. 
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